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SENS, sensitivity analysis Simulation Template With Comments: 

* Instruct the netlist builder to show tolerances \& 0 z 
#tolerance J 

"Suppress automatic vector saves^l ^o* 
tfnosave J 

* Suppress lsSpice4 printout \^ Q ( > 
#noprint S 

"Save vectors needed for measurements! ^ 0 <^ 
#vector J 

"Set the output file pointer to the beginning to remove] 
. * the input net list I 
set rewind -1 

"Set the noecho environment for print formatting ~L 
set noecho J 



"Run the specified simulation and save the results' 

#simulation 

set printmode = save 

#mprint 

"Rename the simulation plotl 
nameplot ref Jt'6 



"Set the print format 

SET COLWIDTH=22 
SET SPICEDIGITS=5 



f)6tOZE 4-/ 



"Tell the user where we are 1 
printstatus -t ########_sensitivity_f or_eachjparameter_########f qzo 



"Loop through all of the parameters 
nextparam null 
while param <> null 

"Alter each parameter 

alterparam tolerance (param) /3 

"Simulate 

#simulation 

Save the parameter reference in the new plot 

paramvec = param 

"Tell the user where we are 

printstatus -p paramvec 

"Save the data 

#mprint 

"Loop through the vectors 
nv = nextvector (null) 
while nv <> null 

*Save the sensitivities of scalar data 

if length (nv) = 1 
nv = nv - ref .nv 
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end; end if 

*Get the next vector 

nv = nextvector (nv) 

end ; end vector loop 

^Restore the parameter 
unalterparam 

*Gef the next parameter 
nextparam 
end; end parameter loop 



*Sef print mode for printing output data 
set printmode = print 
unset noecho 
*Loop through the plots 
pi = nextplot (nuil) 
while pi <> null 

if sameplot (ref .default) = 0 

*Loop through the vectors 
nv = nextvector (null) 
while nv <> null 

if length (nv) = 1 
*Save the sensitivities 
if ref .nv <> 0 

nv = (300*NV) /REF.NV 
else 

nv = 3*NV 

end ; end if 
end ; end if< 
nv = nextvector (nv) 
end ; end vector loop 
end; end if 

*Get next plot 
pi = nextplot (pi) 

end ; end plot loop 

*Tell the user what's going on 

printstatus -t ########sorting_sensitivity_f or_each_j>arameter_## 

printstatus -t 

*Make ref the current plot 

setplot ref 

*Sort by descending value 

sort -vd 

*Loop through the plots 1 

pi = nextplot (null) 
while pi <> null 

if sameplot (ref . default) = 0 

*Print Headers 

SETPARAM PARAMVEC 
printstatus -p paramvec 
ECHO 

ECHO -u ■•**********SENSITIVITY DATA********** 1 
ECHO 

ECHO -un "PARAMETER NAME: 11 
PRINTNAME PARAMVEC 
ECHO 

ECHO -un " NOMINAL VALUE: " 



-4 ZS' 
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PRINTVAL PARAMVEC i F\&uR6 f-3 

ECHO 

ECHO / 
PRINTTEXT -u VECTOR SENSITIVITY% / 
ECHO / 
ECHO 

*Sort by descending data value 
sort -vd 
*Loop through the vectors and print data 
nv = nextvec tor (null) 
while nv <> null 

if length (nv) = 1 
if ref.nv <> 0 

PRINTNAME NV 
PRINTVAL NV 

else 

PRINTNAME NV 
PRINTVAL NV 
ECHO -n * 

end ; end if 
ECHO 

end ; end if 

*Get next vector 

nv = nextvector (nv) 

end; end vector loop 

end ; end if 

*Get next plot 

pi = nextplot(pl) 

end; end plot loop . ________ 

. ECHO •. ' " 

ECHO 

*Print data in output file for SpiceNet tq read _ ? 0 

setplot ref "T- 1 
echo ########## SENSITIVITY analysis Results ############# 
#mprint 
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C:\spice8d\Circuits\Bandpass.cir Setupl — 
*#save V(2) V(3) @C3 [i] @C3 [p] V(l) @Rl[i] @R1 tp] V(4) 
*#save @C4[i] @C4 [p] @R2 [i] @R2 [p] @R5[i] @R5 [p] V(6) V(5) 
*#save@Vl[i] @V1 [p] @V2[i] @V2 [p] @V3 [i] @V3 [p] 
*#alias v_4 v(4n cq± — 
*#view tran v_4 J 
. TRAN ,05ms 20ms 



So? 



50 



C3 
Rl 
C4 
R2 
R5 
XI 



PRINT AC VDB(4) 
OPTIONS vscale=4 
PRINT TRANV4 
ION " 
10K 
ION 
1MEG 
500K 



Sot 



6 5 4 LM324M { } 
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5 
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90 
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DP 
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3 


DX 




EGND 


99 


0 


POLY (2) 3 0 4 0 0 


.5 .5 


FB 


7 


99 


POLY (5) VB VC VE VLP 


VLN 0 53.05E6 


4- -50E6 


50E6 50E6 -50E6 




GA 


6 


0 


11 12 37.70E-6 




GCM 


0 


6 


10 99 11.92E-9 




t? IEE 


3 


10 


DC 2.476E-6 




5 HLIM 


90 


0 


VLIM IK 






11 


2 


13 QX 




< Q2 


12 


1 


14 QX 




y R2 


6 


9 


100. 0E3 




n rci 


4 


11 


26.53E3 




« RC2 


4 


12 


26.53E3 




Z RE1 


13 


10 


4 .820E3 




H RE 2 


14 


10 


4 . 820E3 




R REE 


10 


99 


80.78E6 




ROl 


8 


5 


50 




R02 


7 


99 


50 




& RP 


3 


4 


34 .71E3 




* VB 


9 


0 


DC 0 




VC 


3 


53 


DC 2 




« VE 


54 


4 


DC 5.000E-3 




y VLIM 


7 


8 


DC 0 




*~ VLP 


91 


0 


DC 40 




!: VLN 


0 


92 


DC 40 





Sot 



MODEL DX D(IS=800.0E-18) 
.MODEL QX PNP(IS=800.0E-18 BF=31.58) 
. ENDS 

VI 1 0 AC=1 PULSE 0-1 IMS 
0 DC=-5 
0 DC=5 



V2 5 
V3 6 
.END 
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C : \spice8d\Circuits\Bandpass . Jir Setupl 
. OPTIONS vscale=4 " 
. control 
alias v_4 v(4) 
view tran v^4 
save v(4) foo% 
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set rewind — & )0 
set noecho — £|2- 
TRAN .05ms 20ms "~ 



coo 



set printmode = save 

echo TRAN Analysis Measurements 

echo 

echo Test 1 Mean 
homeCursors 

print Mean (V (4 ) ) 



6>4 



nameplot ref — 6)6 



SET COLWIDTH=22 
SET SPICEDIGITS 



printstatus -t ########_s ensitivity_f or_each_j)arameter__######## 
nextparam null — 
while param <> null 

alterparam tolerance (param) /3 

TRAN .05ms 20ms 



paramvec = param 

printstatus -p paramvec 

echo TRAN Analysis Measurements 

echo 

echo Test 1 Mean 
homeCursors 
print Mean (V (4)) 



5 end 



nv = nextvec tor (null) 
while nv <> null 
if length (nv) = 1 

nv = nv - ref .nv 

end 

nv = nextvector (nv) 
end 

unalterparam 
nextparam 



set printmode 
unset noecho 



print 



pi = nextplot (null) 
while pi <> null 

if sameplot (ref .default) = 0 

nv = nextvector (null) 

while nv <> null 

if length (nv) = 1 
if ref .nv <> 0 

nv = (300*NV) /REF.NV 
else 

nv = 3*NV 

end 

end 

nv = nextvector (nv) 

end 

end 

pi = nextplot (pi) 



end 



printstatus -t ########sorting_sensitivity_f or_each_j?arameter_######## 
printstatus -t £ 2-£> 



setplot ref 
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J 



pl = nextplot (null) 
while pl <> null 

if sameplot (ref .default) = C 

SETPARAM PARAMVEC 

printstatus -p paramvec 



ECHO 

ECHO -u »**********SENSITIVITY DATA**********" 

ECHO 

ECHO -un "PARAMETER NAME: " 

PRINTNAME PARAMVEC 

ECHO 

ECHO -un " NOMINAL VALUE: " 

PRINTVAL PARAMVEC 

ECHO 

ECHO 

PRINTTEXT -u VECTOR SENSITIVITY% 

ECHO 

ECHO 

sort -vd 

nv = nextvector (null) 

while nv <> null 

if length (nv) = 1 
if ref .nv <> 0 

PRINTNAME NV 
PRINTVAL NV 

else 

PRINTNAME NV 
PRINTVAL NV 
ECHO -n * 

end 
ECHO 

end 

nv = nextvector (nv) 

end 



62* 



end 

ECHO 
ECHO 



end 

pl = nextplot (pl) 



setplot ref 

echo ########## SENSITIVITY analysis Results ############# 

echo TRAN Analysis Measurements 

echo 

echo Test 1 Mean 
homeCursors 
print Mean{V (4 ) ) 



.endc 




C3 2 


3 ION 


TOL=5% 


Rl 1 


2 10K 


TOL=2% 


C4 2 


4 ION 


TOL=5% 


R2 2 


0 1MEG TOL=2% 


R5 3 


4 500K TOL=2% 


XI 0 


3 6 5 


4 LM324M { 


.SUBCKT LM324M 12 3 
*■ 


CI 


11 12 


3 .000E-12 


C2 


6 7 


6 . 000E-12 


CEE 


10 99 


315. 8E-15 


DC 


5 53 


DX 


DE 


54 5 


DX 


DLP 


90 91 


DX 


DLN 


92 90 


DX 


DP 


4 3 


DX 


EGND 


99 0 


POLY (2) 3 


FB 


7 99 


POLY (5) VB 
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+ -50E6 
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RP 
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.MODEL DX D(IS=800.0E-18) 

.MODEL QX.PNP(IS=800.0E-18 BF=31.58) 

. ENDS 

VI 1 0 AC=1 PULSE 0-1 IMS 
V2 5 0 DC=-5 
V3 6 0 DC=5 
. END 



******** * *SENS ITI VITY DATA* ********* 



PARAMETER NAME: 
NOMINAL VALUE: 



r5 

500. 00K 



VECTOR SENS ITI VITY% 

mean(v(4)) 1.5111 
******* ***SENSITIVITY DATA********** 



PARAMETER NAME: 
NOMINAL VALUE: 



r2 

l.OOOOMeg 



VECTOR SENSITIVITY% 
mean{v(4) ) 17.265M 
**********SENSITIVITY DATA********** 



PARAMETER NAME: 
NOMINAL VALUE: 



C4 

10.0000N 



VECTOR SENSITIVITY% 
mean(v(4) ) -752. 77M 

****** ****SENSITIVITY DATA********** 



PARAMETER NAME: 



rl 

10.0000K 



NOMINAL VALUE: 

VECTOR SENSITIVITY% 
mean(v(4)) -571. 46M 

**********SENSITIVITY DATA********** ✓ 



PARAMETER NAME: 
NOMINAL VALUE: 



VECTOR 
mean (v (4) ) 



c3 

10.0000N 

SENSITIVITY% 
4.5201 
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Fl6rORE 7 



70 O 



########## sensitivity analysis results ############# 
tran analysis measurements 



\4 test 1 mean 

mean(v(4)) = 2 . 086052e-001 

Total run time: 0.583 seconds. 

Total run time: 0.583 seconds. 

Memory remaining = 1996210 Kbytes 
Memory Used = 14401 Kbytes 



13 of 24 




1 5977-1 3 



RSS, root summed square analysis Simulation Template With Comments: 



*lnstruct the netlist builder to show tolerancesl <* o2 
#tolerance J 



* Suppress automatic vector savesl 
#nosave J 



* Suppress lsSpice4 printout \ <% 0 £ 
#noprint J 

*Save vectors needed for measurements! ?o g 

#vector J 

*Set the output file pointer to the beginning to remove] 

* the input net list ' fi ° 
set rewind 




*Set the noecho environment for print formatting! g ^ 
set noecho J 



*Run the specified simulation and save the results 

#simulation 

set printmode = save , 

#mprint 



*Set the print format 

SET COLWIDTH=22 
SET SPICEDIGITS=5 



*Rename the simulation plot 
nameplot ref 



876 



*Loop through all of the parameters 
nextparam null « 

*Tell the user where we are ] 

printstatus -t »######## sensitivity for -each parameter ######## 
while param <> null ( 

* Alter each parameter * 
alterparam tolerance (param) /3 
"Simulate 
#simulation 

*Save the parameter reference in the new plot 
paramvec = param 

* Tell the user where we are I — % 22 
printstatus -p paramvec 
*Save the data 
#mprint 

"Loop through the vectors 
nv = nextvector (null) 
while nv < > null 
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"Save the sensitivities of scalar data 
if length (nv) = 1 
nv = nv - ref.nv 

end ; end if 
"Get the next vector 
nv = nextvector (nv) 
end ; end vector loop 
"Restore the parameter 
unalterparam 
*Get the next parameter 
nextparam 
end ; end parameter loop 
*Set print mode for printing output data 
set printmode = print 
unset noecho 

*Make a new plot for rss data and make it the current plot 
newplot rss ref. default ref. default 
setplot rss 
*Loop through the vectors 
nv = nextvector (null) 
while nv <> null 

"Initialize the scalar data to zero 

IF LENGTH (NV) = 1 
nv = 0 

end ; end if 

" Get the next vector 
nv .= nextvector (nv) 

end ; end vector loop 

*Loop through the plots 
pi = nextplot (null) 
*Tell the user where we are 

printstatus -t »•######## computing rss for each measurement #### 
while pi <> null 
*Select plots 

if sameplot (rss .default) = 0 

if sameplot (ref .default) = 0 



S32 



9Zb 



"Tell the user what we are doing 
printstatus -1 pi 

*Make the saved parameter reference, paramvec current 

SETPARAM PARAMVEC 

*Gef the next vector 

nv = nextvector (null) 

"Print formatted data 

ECHO 

ECHO -U ''**********RSS DATA**********" 
ECHO 

ECHO -UN " PARAMETER NAME: » 

PRINTNAME PARAMVEC 

ECHO 

ECHO -UN " NOMINAL VALUE: » 

PRINTVAL PARAMVEC 

ECHO 

ECHO -UN "TOLERANCE VALUE: " 
PRINTTOL PARAMVEC 
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ECHO _ ^ 0 

ECHO YlZon 
PRINTTEXT -UN VECTOR 

PRINTTEXT -U SENSITIVITY% RSS_CONTRIBUTION 
ECHO 
ECHO 

*Loop through vectors 
while nv <> null 

if length (nv) = 1 
IF REF .NV <> 0 

*Calculate the RSS percentage if value is not zero 

PRINTNAME NV 
NEWNV = (300*NV) /REF.NV 
PRINTVAL NEWNV 

ELSE 

*Calculate the RSS if value is not zero 

PRINTNAME NV 
NEWNV = 3*NV 
PRINTVAL NEWNV 
ECHO -N * 

end ;end if 

*Save and print the sum of squares 
rss.nv = rss.nv + nv * nv 
PRINTVAL RSS.NV 
ECHO 
end ;end if 
* Get next vector 
nv = nextvector (nv) 
end ;end vector loop 
* Sort plot by descending value 
sort -vd 
end ;end if 
end ;end if 
* Get next plot 

pi = nextplot(pl) 

end ;end plot loop 

*Sort the rss plot by descending value 
setplot rss 
SORT -VD 
*Print Headers 

ECHO 

ECHO -U ********** *RSS HI/LO ANALYSIS RESULTS**********" 
ECHO 

PRINTTEXT -UN VECTOR 

SET COLWIDTH=15 ' K # 

PRINTTEXT -U NOMINAL RSS-VALUE TOLERANCE% HI_VALUE LO__VALUE 
ECHO 
ECHO 

*Ma/ce a new plot for results 

newplot hirss ref. default ref. default 

*Loop through the vectors 
nv = nextvector (null) 
while nv <> null 

if length (nv) = 1 
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*Print formatted data 



SET COLWIDTH=22 

PRINTNAME NV 
SET C0LWIDTH=15 
PRINTVAL REF .NV 

NV = 3 * SQRT (ABS (NV) ) 
PRINTVAL NV 

IF REF.NV <> 0 

NEWNV1 = ( 10 0*NV) /REF.NV 

ELSE 

NEWNV1 = NV*0 

END 

PRINTVAL NEWNV1 

NV = REF.NV + NV 
HI_RSS = REF.NV + NV 
LO_RSS = REF.NV - NV 
PRINTVAL HI_RSS 
PRINTVAL LO_RSS 
ECHO 

end ; end if 

*Get next vector 

nv = nextvector (nv) 



end; end vector loop t 

ECHO _ 
ECHO 

*Print data in output file for SpiceNet to read 
setplot hirss 

echo ########## RSS HI analysis Results ############# 
ttmprint 

RUSAGE ELAPSED 
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EVA, Extreme Value Analysis Simulation Template With Comments: 

* Instruct the netlist builder to show tolerances 1 
#tolerance J 

* Suppress automatic vector saves 1 



ttnosave 

"Suppress lsSpice4 printout 
ttnoprint 



1 



"Save vectors needed for measurements 
#vector i. 



Set the noecho environment for print formatting 
set rewind — ^10 
set noecho — 

"Run the specified simulation and save the results 
#simulation 

pltype = 0 Jdentify the plot type for later use r °\ 

set printmode = save 
#mprint 

"Set the print format 

SET COLWIDTH=22 - 
SET SPICEDIGITS=5J 
nameplot ref - q/£r 

newplot evahi ref -default ref. default 1 
evahi . pltype = 0 Jdentify the plot type for later use 
"Print status for the user 

printstatus -t "######## sensitivity for each parameter ########" 

"Loop through the parameters 

nextparam null 
while param <> null 

"Alter each parameter 

alterparam tolerance (param) /3 

"Simulate, making a new plot for results 

#simulation 

"Save the current parameter reference 
paramvec = param 
*7e// the user where we are 
printstatus -p paramvec 
"Save the tol and paramval 
paramtol = tolerance (param) 
paramval = getparam (param) 
pltype = 1 Jdentify the plot type for later use 
"Save the simulation results 
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#mprint 

*Loop through all the vectors 
nv = nextvector (null) 
while nv <> null 

"Save the sensitivities for all scalar measurements except pltype 
if length (nv) = 1 

if nv <> pltype 

nv = nv - ref.nv 
end ;end if 
end ;end if 

nv = nextvector (nv) 

end ;end vector loop 
*Restore the param 
unaiterparam 
*Get the next param 
nextparam 
end;end parameter loop 



fi6ot?£ 1-2 



H . a 



rn 



*Make ref the current plot 

setplot ref 

*Loop through the vectors in ref 

nv = nextvector (null) 
"Tell the user where we are 

printstatus -t "######## measurements ########" 
while nv <> null , 
if length (nv) = 1 

*Loop through all the plots containing scalar vectors 
pi = nextplot (null) 
while pi <> null 
if pltype = 1 

* the inner loop, we are looping through each sensitivity plot looking at the same vector 

* we will alter the parameter id'd byparamvec to maximize/minimize the vector 

setparam paramvec 

"Change each parameter to its worst case extreme value 
if nv >= 0 

alterparam paramtol 

else 

alterparam -paramtol 

end 

end 

pi = nextplot (pi) 

end 

"Simulate for the extreme case and save the data in a new plot 

#simulation 
#mprint 

pltype = 2 ^Identify the plot type for later use 

* if we want sensitivity at the extreme, we weed to go through each param 

* and change it to be a bit different than it is at the extreme, run a simulation , 

* and mark it as pltype 3 along with its paramvec, then we can take the diff 

* from the pltype = 2 to get the sensitivity at the extreme if the sign at the extreme 



19 of 24 



□ 



15 



)77-13 



* /s different than at the nominal, we can report an error or go on to do worst case 

* for worst case, we need to reduce the param change by 1/2 and do this over again... 

* either continue in this loop or make a wc loop afterward... save the paramvalue 

* and tolerance 

evahi . nv = nv 
if nv <> pi type 

*Tell the user where we are 
printstatus -n nv 
end ;endif 

end ;end plot loop 
iget the next vector 

nv = nextvector (nv) 

end ;end vector loop _____ - » _ 

*Set print mode and print header — 

set printmode = print 
unset noecho 
setplot EVAHI 
ECHO 

ECHO -U ''**********EVA PARAMETER LIST**********" 
ECHO 

PRINTTEXT -UN PARAMETER 
PRINTTEXT -U NOMINAL TOLERANCE 
ECHO 
ECHO 

*Loop through the parameters 
nextparam null 
while param <> null 

^Extract the saved param reference and print its data 

paramvec = param 

PRINTNAME PARAMVEC 

PRINTVAL PARAMVEC 

PRINTTOL PARAMVEC 

ECHO 

*get the next param 
nextparam 

end ;end param loop 

echo ;print a blank line — 

*Make a new plot to hold sorted results 
newplot evasort ref. default ref. default 
*Make ref the current plot 
setplot REF 

*Loop through all vectors in ref 
nv = nextvector (null) 
while nv <> null 

if length (nv) = 1 

*save the result in evasort as a percent of its value 
if ref.nv <> 0 

evasort .nv = ( (evahi . nv-ref . nv) *100 ) /ref . nv 
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else 

evasort.nv = 0; 

end ;endif 

end ;end if 
*Get the next vector 

nv = nextvector (nv) 

end ;end vector loop 
*Print some headers 

ECHO 

ECHO -U "****** ****EVA-HI RESULTS 
ECHO 

PRINTTEXT -UN VECTOR 

PRINTTEXT -U NOMINAL EVA-HI CHANGE % 

ECHO 

ECHO 

setplot evasort 

*Sort evasort by descanting data 

sort -VD 

"Loop through the vectors 
nv = nextvector (null) 
while nv <> null 

*lf its the correct data in the correct plot, print it 
if length (nv) = 1 

if nv <> pltype , 
PRINTNAME NV 
PRINTVAL 
PRINTVAL 
PRINTVAL 
ECHO 

end ; end if 

end ; end if 

*Get the next vector 

nv = nextvector (nv) 
end ; end vector loop 
ECHO 
ECHO 



* now the eva results are in pltype - 2 plots 
*Print the results so SpiceNet can read the eva-hi data 

set printmode = print 
unset noecho 
setplot evahi 

echo ########## EVA HI analysis Results ########.##### 
#mprint 



REF.NV 
EVAHI.NV 
EVASORT . NV 
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WCS, Worst Case by Sensitivity Simulation Template With Comments: 



* Instruct the netlist builder to show tolerances"} 
#tolerance 



r 



"Suppress automatic vector saves! 



#nosave 



* "Suppress lsSpice4 phntout'T )n ^, 
#noprint J IUOh 

"Save vectors needed for measurements! 




#vector 



\00*€ 



"Set the output file pointer to the beginning to remove] 

" the input net list J 1010 

set rewind 

"Set the noecho environment for print formatting! j0iz 
set noecho J 



*Run the specified simulation and save the results 
#simulation 



set printmode = save 
#mprint 

"Set the print format 

SET COLWIDTH=22 f- \6\% 
SET SPICEDIGITS=5 



- 1014 



*Rename the simulation plotL 
nameplot ref J 



1016 



"Make a newplot for results I 
newplot result ref .default ref. default - 

"Set the current plot to ref 1 
setplot ref J w 

"Print status for the user 

printstatus -t »######## sensitivity for each parameter ######## L 



"Loop through the parameters 
nextparam null 
while param <> null 
"Alter each parameter 
alterparam tolerance (param) /3 
"Simulate, making a new plot for results 

#simulation 
"Save the current parameter reference 

paramvec = param 
"Inform the user about what's being done 

printstatus -p paramvec 
"Make and save the measurements 
#mprint 
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*Loop through the vectors 

nv = next vector (null) ' 
while nv <> null 
*Save the sensitivity of scalar quantities 
if length (nv) = 1 
nv = nv - ref.nv 
*Save the worst case -hi value 
result. nv = result. nv + abs(3*nv) 
end \endif 

nv = nextvector (nv) 
end \end vector loop 
* restore the parameter value 
unalterparam 
*gef the next parameter 
nextparam 

end -end parameter loo p 

*Set the print mode to print instead of save 

set printmode = print 

* Restore the echo mode for printing 

unset noecho 

*Set result to the current plot 

setplot result 

* Print the header 

ECHO 

ECHO -U »»**********WCS PARAMETER LIST**********" 
ECHO 

PRINTTEXT -UN PARAMETER 
PRINTTEXT -U NOMINAL TOLERANCE 
ECHO 
ECHO 

*Loop through the parameters 
nextparam null 
while param <> null 

paramvec = param 



ECHO ' 
nextparam 

end 

ECHO J 

*Make a new plot to hold sorted results 

newplot wcsort ref. default ref . default 

*Set the current plot to ref 
setplot REF 

*Loop through its vectors 
nv = nextvector (null) 
while nv <> null 

* Calculate the wc as a percent change results 
if length (nv) = 1 
if ref.nv <> 0 

wcsort. nv = ( (result . nv-ref . nv) *100 ) /ref . nv 



* Print the row 



PRINTNAME PARAMVEC 
PRINTVAL PARAMVEC 
PRINTTOL PARAMVEC 
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else 

wcsort.nv = 0; 
end ;end if 
end ;end if 

nv = nextvector (nv) 

end ; end vector loop 



\ 
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"Print headers 

ECHO 

ECHO -U »**********WCS-HI RESULTS**********" 
ECHO 

PRINTTEXT -UN VECTOR 

PRINTTEXT -U NOMINAL WCS-HI CHANGE % 
ECHO 
ECHO 

;sort wcsort by descending value 
setplot wcsort 
'sort -VD 
"Print the ordered list 
nv = nextvector (null) 
while nv <> null 

if length (nv) = 1 

PRINTNAME NV 

PRINTVAL REF.NV 

PRINTVAL RESULT . NV 

PRINTVAL WCSORT . NV 
ECHO 

end 

nv = nextvector (nv) 

end 
ECHO 

ECHO _ " " : 

*Set the current plot to the wc results . 
setplot result 

echo ########## WCS HI analysis Results ############# 
"Print the measured results in a form that can be read back 
"into SpiceNet 
#mprint 

"Report the elapsed time in the output file 

rusage elapsed — 



- /0 3Q 
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